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LINEAR  THEORY  OF  THE  RAYLEIGH-TAYLOR  INSTABILITY  AT  THE 
INTERFACE  BETWEEN  TWO  COMPRESSIBLE  MEDIA 

1.  INTRODUCTION' 

The  Rayleigh-Taylor  instability  occurs  when  a  fluid  supports  a  denser 
fliid  against  gravity,  whereupon  the  two  tend  to  interchange  positions.  It 
is  encountered  frequently  in  nature  and  in  the  laboratory.  In  inertial 
confinement  fusion  experiments ,  in  which  an  ablatively  driven  median 
implodes,  compressing  the  material  ahead  of  the  ablation  front  to  high 
densities,  Rayleigh-Taylor  instabilities  can  play  a  crucial  role.  They  may 
occur  in  the  ablative  region,  at  the  compression  front,  or  (in  the  case  of 
layered  targets)  at  an  interface  between  layers  of  different  density. 

The  familiar  treatments  of  the  Rayleigh-Taylor  instability  which  give 
rise  to  the  classical  dispersion  relation  (growth  rate  proportional  to  the 
square  root  of  the  product  of  perturbation  wavenumber  k  and  gravitational 
acceleration  g)  presuppose  that  the  fliid  moves  incompressibly ,  i.e. ,  that 
the  iivergence  of  the  flow  field  vanishes.  It  is  frequently  argued  that  when 
the  effects  of  material  compressibility  are  taken  into  account,  the  growth 
rates  will  re  reduced  in  comparison  with  the  incompressible  case,  since 
compression  absorbs  some  of  the  energy  which  would  otherwise  go  into  fluid 
motion.  One  might  speculate,  for  example,  that  growth  rates  at  long 
wavelengths  would  be  limited  to  a  quantity  of  order  the  speed  of  sound 
divided  by  the  wavelength.  However,  Bernstein  et  al. 1  have  shown  that  in  a 
broad  :iass  of  general  hydromagnetic  systems  the  unstable  moles  with  lowest 
threshold  are  associated  with  incompressible  perturbations.  It  is  thus 
conceivable  that  a  compressible  system  might  exhibit  incompressible  modes  of 
interchange  instability  with  the  classical  Rayleigh-Taylor  growth  rate. 
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This  is  in  fact  the  case  for  it  least  seme  models  of  imploding  hy droiiynariic 
systems.-*  *  To  late  r.o  one  seems  to  have  suggested  that  compressibility 
might  enhance  the  instability  growth  rates. 

Previous  attempts  to  determine  the  effect  of  compressibility  on  the 
Payleigh -Taylor  instability  have  tended  to  include  numerous  other  effects  as 
well,**  obscuring  the  role  played  by  compressibility ,  or  have  neglected  to 
state  clearly  the  assumptions  used  in  the  derivation,^’0  leaving  in  doubt 
its  domain  of  validity.  The  only  conclusicr.  these  authors  appear  to  have 
reached  regarding  compressibility  is  that  its  effects  are  small.  In 
contrast,  ve  have  chosen  to  treat  a  model  err  clem  w.nich  can  be  solved 
exactly.  We  analyse  the  simple  case  of  two  contiguous  se.mi -infinite  slab 
regions  in  a  constant  gravitational  field  each  filled  with  an  ideal  fluid 
whose  density  decreases  exponentially  in  the  upward  direction  (Fig.  1).  If 
we  subject  this  state  to  an  infinitesimal  certurbat ion,  the  resulting  problem 
i  i  free  of  such  complicating  features  as  thermal  conduction,  time  depenler.ce 
in  the  unperturbed  state,  curvilinear  geometry,  material  flew  across  or  along 
the  interface,  etc.,  and  can  ce  solved  rigorously  for  the  eigenfre  juer.cy  v 
for  either  compressible  (pclytropic)  or  incompressible  fluids.. 

When  the  media  are  incompressible  ve  recover  in  the  short-wavelength  limit 
the  expected  result  uj2  =  -kg(  P0-Pg' )/(  P,.+  P3' ) ,  where  Pq  is  the  unperturbed 
density  of  the  upper  medium  and  p^’  is  that  of  the  lover,  both  measured  an 
the  interface.  For  longer  wavelengths  this  simple  form  is  replaced  by  one  in 
which  the  respective  scale  heights  p2/pcg  and  p0/p0'g  enter,  where  pc  is  the 
unperturbed  pressure  at  the  interface.  Ir.  the  limit  k  +  0,  u>  becomes 
proportional  to  k.  This  wavenumber  dependence  reflects  the  structure  of  the 
basic  state.  In  a  (possibly  more  realistic)  model  involving  two  layers  of 
finite  thickness  d,  d' ,  the  dispersion  relation  would  n.aiural1y  depend  or.  the 
pientities  kd,  hi'  and  would  be  still  more  complicated. 
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Following  the  treatment  of  the  incompressible  case  we  tarn  to  she 


comcresa'!  cle  case. 


the  short -wavelength  limit  we  cnee  again  obtain  the 


classical  result.  This  is  unsurprising,  as  can  depend  on  the  sound 

speeds  c,  c'  only  through  the  ; i mens i on les s  ratios  kc2/g,  kc'^/g.  barge 
values  of  k  are  thus  equivalent  to  large  values  of  the  sound  speed,  i.e.  ,  to 
incompressibiltv .  When  .<  is  finite,  however,  we  unexpectedly  find  that  the 
growth  rates  are  larger  than  the  incompressible  ones.  We  illustrate  the 
result  by  solving  the  dispersion  relation  numerically  as  a  function  of 


wavenumber  for  two  different  values  of  the  density  jump  at  the  interface 


for  several  different  '/a lues  of  y,  assuming  that  the  latter  has  the  same 
value  in  both  regions.  The  difference  in  u)2/kg  as  y  varies  between  unity  an 
00  car.  amount  to  a  factor  -  2.  These  results,  though  they  do  not  constitute 
an  exhaustive  parameter  survey,  typify  what  may  be  e-.pected  in  real  physical 
situations. 

We  conclude  the  paper  with  a  brief  discussion  of  cur  findings  ar.d  seme 
speculations  on  their  significance. 
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INCOMPRESSIBLE  FLUIDS 


'cnsider  a  fluid  described  by  the  continuity  equation 

_l£.  +  7 .  ov  +  D 
3t  ~ 


(1) 


a n trie  momentum  equation 


P  (-ry  +  v»7v)  =  -7p  -  pge^. 


(2) 


/here  g  is  a  constant.  Thus  -gev  is  a  uni  fern  gravitational 


acceleration,  oriente^ 


: he  downward  (ns 


/ '  direction.  These 


equations  have  an  equilibrium  solution 

v  =  0; 

<  _  -> 
p/p  =  c- 


const ; 

7n  =  -pge  =  -e  c2-^-, 

"ty  -y  3y 


'3) 

(U) 

f c  'i 


-l.r\  -  n  _v-  /  •>  2  \  '<C> 

p..y  ,  -  P0-Xt.l -bj  /-  >  >  'c- 

P3  constant.  We  take  the  interface  separating  the  two  media  in  the  proble 
to  be  located  at  y  =  C,  i.e.  ,  coinciding  with  the  x-axis  (see  rig.  l).  To 
distinguish  the  two  regions  we  'will  label  all  quantities  belonging  to  the 
lower  one  with  primes.  Mechanical  equilibrium  requires  pressure  balance, 

p(0)  =  p* (O),  ( 7 ) 

but  in  general  pq  *  pQ’  and  c  4  c'.  Since  we  anticipate  that  instability 
will  occur  when  the  upper  medi’un  is  denser  (at  y  =  C)  than  the  lower,  we 
assume  o'  >  c. 

Suppose  that  this  state  is  subjected  to  an  inf initesimal  perturbation 
defined  by  the  local  displacement  £(x,y,t)  of  an  element  of  fluid. 

Using  the  subscript  "1"  to  distinguish  perturbed  quantities,  we  have  for  t 
velocity 
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,,  „2  35., 

3  A  Z  J 

+  z  — —  -  — 2  ~ — 

~  dv  C  c  di 


..et  ,13  3ee.v 


c:  the  :or!r* 


cp(  iux  -  10'  -  lot), 


"  expU-ux  -  qy 


equations  (2C;-(21)  include  a  factor  exp  (-i«rc),  Walsh  nay  be  either 
oscillatory  or  exponential,  depending  on  v.  Its  presence  implies  tr.at  vs  ar  = 
meeting  ordinary  (nonsingular J  .aortal  modes  of  the  system.  Unis  restriction 
is  equivalent  to  solving  the  initial  value  problem  associated  with  -q.  (i-0, 
e.g. ,  by  carrying  cut  a  Laplace  transform,  and  then  throwing  out  the  portion 
of  the  solution  which  depends  on  the  initial  data.  lie  return  to  this  point 
in  our  concluding  discussion  (Section  u). 

Substitution  of  Tqs.  (20)-(2l)  in  Eqs.  (l3)-(l9)  yte_is 

-<o2A  =  h=23  +  g3  -  (g2/c2\A,  (22) 

0  =  -c2(u2-h2)B  -  gu3  +  (g2/c2)uA.  (2c) 

Thus  u  can  be  found  from  the  solubility  condition,  i.e. ,  by  eliminating  A/3 


A  _  gu  -  c  2(k2-u2 ) 

3  ug^/c2 


A  uc 2  + 


.vi.ng  the  resulting  quadratic  equation,  ve  find  in  the  upper  half  plane. 


taxing  the  square  root  with  positive  real  part 
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Fig.  2  —  Dimensionless  squared  growth  rate  —  <*.>-  log  vs  wavpnumber  k 
for  the  two  basic  states  shown  in  Fig.  1(a),  (b).  with  the  same  choice 
of  units  as  in  the  latter.  The  adiabatic  index  y  in  both  regions  is  taken 
to  be  1,  5/3  or  «=,  as  indicated  by  the  label.  Note  that  the  curves  as¬ 
ymptotically  approach  the  value  >p  —  o')  ip  +  p ' ),  equal  to  0.333  and 
0.818.  respectively. 


U) 


CONCLUSION'S 
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F.vi  ient’v  the  onset  of  the  process  by  which  the  two  participating  fl 
Interchange  positions  is  ailed  by  the  extra  degrees  of  freedom,  which  alio 
then:  to  compress  and  expand.  A  mechanical  analogy  is  supplied  by  the  e/.a 
of  the  inverted  physical  pendulum  (e.g.  ,  a  brick  smokestack ) ,  which  falls 
over  upon  undergoing  on  angular  displacement.  different  parts  of  the 
pendulum  falling  individually  would  assume  different  angular  velocities  u 
to  the  action  of  gravity  alone.  Hence  the  pendulum  tends  to  bend  or  brea 
under  internal  stresses,  allowing  it  to  fall  over  faster  than  or.e  constra 
to  remain  rigid. 

Since  the  range  of  effective  values  of  y  for  materials  of  interest  i 
inertial  confinement  fusion  is  not  great,  the  quantitative  changes  in 
instability  growth  rates  (vis-a-vis  those  found  in  incompressible  theory) 
like  ly  to  be  modest,  perhaps  a  few  tens  of  percent.  7urth.erm.ore,  the  pre 
calculation  is  admittedly  idealised,  and  the  results  are  certainly  model- 
dependent.  Nevertheless ,  our  findings  do  not  encourage  or.e  to  anticipate 
that  compressibility  effects  can  mitigate  the  Hay le i gh -Taylor  instability 
Perhaps  the  only  conclusion  should  be  that  the  effects  of  compressibility 
need  to  be  carefully  treated  in  each  individual  situation. 
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